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CRITICAL COMBINATIONS OF SHEAR AND TRANSVERSE DIRECT STRESS FOR AN INFINITELY 
LONG FLAT PLATE WITH EDGES ELASTICALLY RESTRAINED AGAINST ROTATION 

By S. B. Ba.tx>orf and John C. Hotjboht 


SUMMARY 

An exact solution and a closely concurring approximate 
energy solution are given for the buckling of an infinitely long 
flat plate under combined shear and transverse direct stress 
with edges elastically restrained against rotation. It was found 
that an appreciable fraction of the critical stress in pure shear 
may be applied to the plate without any reduction in the trans- 
verse compressive stress necessary to produce buckling. An 
interaction formula in general use was shown to be decidedly 
conservative for the range in which it is supposed to apply. 

INTRODUCTION 

In the design of stressed-skin structures, consideration 
must sometimes be given to tho critical stresses for a sheet 
under a combination of shear and direct stress. The upper 
surface of a wing in normal flight, for example, is subjected 
to combined shear and compressive stress, and the lower 
surface is subjected to combined shear and tensile stress. 
The upper surface may then buckle at a lower compressive 
or shear stress than if either stress were acting alone. The 
critical shear stress for the lower surface may be increased 
by the presence of the tensile stress. 

If the wing has closely spaced chordwise stiffeners, the 
skin between two adjacent stiffeners may be regarded as a 
long sheet slightly curved in the longitudinal or chordwise 
direction, straight in the transverse or spanwise direction, 
and loaded in shear and transverse direct stress. A con- 
servative preliminary estimate of the critical stresses may be 
obtained if the sheet is considered to be flat and infinitely 
long. In the present paper, the critical stresses are computed 
for an infinitely long flat plate loaded as indicated in figure 
1 (a). The corresponding idealization of the case of a wing 
with spanwise stiffeners, show in figure 1 (b), was treated 
in reference 1. 

CONVENTIONAL INTERACTION FORMULAS 

For buckling of structures under combined loading con- 
ditions, no general theory has been developed that is appli- 
cable to all cases. Stress ratios (reference 2), however, 
provide a convenient method of representing such conditions. 
For example, the ratio of the shear stress actually present 
in a structure to the critical shear stress of the structure 
when no other stresses are present may be called the shear- 
stress ratio. Stress ratios may similarly be defined for each 
type of stress occurring in the structure. 



(b) 


fa) Type of loading problem solved In present paper. 

(b) Type of loading problem solved In reference 1. 

Figure 1. — Buckling of an Infinitely long plate nnder combined loads. 

It is generally assumed that equations of the type 

... =1 ( 1 ) 

may be used to express the buckling conditions in the case 
of combined loading (reference 3, pp. 1-18) . In equation (1) , 
i?i, Ri, and R 3 are stress ratios and p, g, and r are exponents 
chosen to fit the known results. (All symbols are defined in 
appendix A.) Such a formula gives the correct results when 
only one type of loading is present and has the further 
advantage of being nondimensional. Equation (1) implies, 
moreover, that the presence of any positive fraction of the 
critical stress of one type reduces the amount of another 
type of stress required to produce buckles; this implication 
appears reasonable and has been proved time in some cases 
(references 1, 2, and 4). 

In reference 2 the following interaction formula is given 
for an infinitely long plate with clamped edges loaded in 
shear and longitudinal compression: 

Rf s +Rc= 1 ( 2 ) 

where R, is shear-stress ratio and R e is longitudinal direct- 
stress ratio. The same formula is recommended in reference 
3 for general use for the buckling of any fiat rectangular 
plate, regardless of the direction of compression and the 
degree of edge restraint. 

Later theoretical work (reference 1) shows that, to a high 
degree of accuracy, for an infinitely long plate with any 
degree of edge restraint loaded in shear and longitudinal 
compression 

R, 2 +R e = 1 (3) 
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Figure 2.— Crltloal combinations oi shear-stress coefficient k, and transverse direct-stress coefficient A. for an Infinitely long plate with edges elastically restrained against rotation. 
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The same formula was found in reference 5 to be applicable 
to simply supported rectangular plates of aspect ratios 

0.5, 1, and 2; the conclusion was drawn that interaction 
curves in stress-ratio form are practically independent of 
the dimensions of the plate. 

The present analysis, however, indicates that the buckling 
of an infinitely long plate loaded in shear and transverse 
compression (fig. 1 (a)) is not adequately represented either 
by equation (2) or (3) or by any formula of the type of 
equation (1). Two independent theoretical solutions to 
this buckling problem are given in appendixes B and C. 
Appendix B contains the exact solution of the differential 
equation of equilibrium, and appendix C contains an energy 
solution leading directly to an interaction formula. This 
energy solution, which gives approximate values only, was 
made to obtain an initial quick survey of the problem and 
to provide a check on the results of the exact solution. 
Approximate interaction formulas in substantial agreement 
with these results were given for the cases of simply sup- 
ported and clamped edges in reference 6. 

RESULTS AND DISCUSSION 

In figure 2, curves are given that indicate the critical 
combinations of shear and transverse direct stress for an 
infinitely long plate with edges elastically restrained against 
rotation. These curves are computed from the exact solu- 
tion presented in appendix B. The degree of edge restraint 
is denoted by e, which is defined in appendix B in such a 
way that zero edge restraint corresponds to simply supported 
edges and infinite edge restraint indicates clamped edges. 
A similar set of curves is given in terms of stress ratios in 
figure 3. The numerical values used to plot-figures 2 and 3, 
together with the values found by the energy solution, are 
given in table I. 

The most striking feature of these results is that an appre- 
ciable fraction of the critical stress in pure shear can evi- 


dently be applied to the plate without any reduction in the 
compressive stress necessary to produce buckling. (See 
fig. 3.) This fraction varies from about one-third to more 
than one-half, depending on the degree of restraint. At 
shear stresses higher than those corresponding to this frac- 
tion, the compressive stress required to produce buckling is 
reduced by the presence of shear. The result that the com- 
pressive buckhng stress is entirely unaffected by the presence 
of a considerable amount of shear is probably peculiar to 
infinitely long plates. It is to be expected, however, that 
this result will be closely approached in the caso of long 
finite plates. 

In figure 4, a comparison is made between the exact solu- 
tions und the interaction formulas of equations (2) and (3). 
Equation (2), which is the interaction formula in general 
use, is seen to be decidedly conservative. 

CONCLUSIONS 

The exact solution of the differential equation for the 
buckling of an infinitely long flat plate under combined shear 
and transverse direct stress with edges elastically restrained 
against rotation indicates the following: 

1. An infinitely long flat~plate may he loaded with an 
appreciable fraction of its critical stress in pure shear without 
causing any reduction in the transverse compressive stress 
necessary to produce buckling. 

2. An interaction formula in general use for rectangular 
plates in combined shear and compression is decidedly con- 
servative when applied to an infinitely long plate in shear 
and transverse compression. 


Langley Memorial Aeronautical -Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., November 8, 1944~ 


















APPENDIX A 

SYMBOLS 


Ci, Cs functions of edge restraint coefficient e given in 
appendix C 

D flexural stiffness of plate per unit length, in-lb 



elastic modulus of material, psi 
compressive force per unit length, lb/in. 
shearing force per unit length of plate, Ib/in. 
rotational stiffness per inch of restraining member 
at edge of plate, lb/radian 

T c work done by compressive force per half wave 

length, in-lb 

T, work done by shear force per half wave length, 

in-lb 

Vi strain energy in plate per. half wave length, in-lb 

V 2 strain energy in edge restraint per half wave length, 

in-lb 

7 function of y associated with deflection of plate 

during buckling 

b width of plate, in. 

bi width of plate in oblique coordinate system of ref- 

erence 1, in. 

c function of a, j 3 , and X 

k c , k, critical compressive and shear-stress coefficients, 
respectively 


m root of a characteristic equation of appendix B 

y,,y 2 ,/s functions of restraint coefficient e given in 
appendix C 

t thickness of plate, in. 

wb 

U= 2X 

w displacement of buckled plate from original position 

w 0 amplitude of assumed wave form of buckle 

x longitudinal coordinate of plate 

y transverse coordinate of plate 

a, p functions of X, 7, and k c 

e nondimensional coefficient of edge restraint 

7 one of two parameters determining buckle form 

X one of two parameters determining buckle form 

(half wave length of buckle, in.) 

[x Poisson’s ratio 

a direct stress 

<t v transverse direct stress, psi 

r shear stress, psi 

<j> angle between buckle node and y-axis 

0=tan 0 

Ri, R t , Rg stress ratios 
R, shear-stress ratio 

R e longitudinal direct-stress ratio 

cr critical (used as subscript) 
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APPENDIX B 

SOLUTION BY DIFFERENTIAL EQUATION 


Statement of problem. — The exact solution for the critical 
stress at which buckling occurs in a flat rectangular plate 
subjected to combined shear and compression in its own 
plane may be obtained by solving the differential equation 
that expresses the equilibrium of the buckled plate. The 
plate is assumed to be infinitely long, and equal elastic 
restraints against rotation are assumed to be present along 
the two edges of the plate. 

Differential equation. — Figure 5 shows the coordinate 
system used. The differential equation for equilibrium of a 
flat plate under shear and transverse direct stress is (from 
reference 7) 

w-° <=» 

It is convenient to write <r„ and r in terms of the dimension- 
less buckling coefficients k„ and k, by means of the relations 


k c Tr 2 D \ 

= bH 

k.^D 

= m 


(B2) 


Substitution of the expressions for a, and r from equations 
(B2) in equation (Bl) gives 

d*w , „ d*w , d*w , 2tc*k, tfw , <¥w_ n . 

da: i ^ 2 dx J dj/ 4 ~^dy J ' b* drdy b* dy* x ' 

Solution of differential equation. — If the plate is infinitely 
long in the ^-direction, all displacements must be periodic 
in x and the deflection surface may be taken in the form 


£2 

w=7e x 


(B4) 


where Fis a function of y only and X is the half wave length 
of the buckles in the r-direction. 


< 

> 

|— ~ t — .r-r — ^ 

b 

T 

* 

1 

j 


l ] 



rrrrnrrr 


t-T-'l-f 


Fracax 5,—C o ordinate system used fn appendix B for an Infinitely long plate under combined 
shear and transverse compression. 


Substitution of the expression for w from equation (B4) 
in the differential equation (B3) gives the following as the 
equation that determines 7: 

<77 , /x 2 £ c 2v*\d?7 , 2x 3 ik,d7 , x* xr „ /r>lrS 

"° (B6) 

A solution of equation (B5) is 

TT 

7—e b 

where m is a root of the characteristic equation 

m<+[ 2 (yJ“ ***•] w*~ 2 x 3 (y) k,m+(~y=0 (B6) 

Except for the substitution of 2 ~t?k a for 2 ’ 

equation (B6) is identical with equation (A-6) of appendix 
A of reference 8, in which equation (B3) of this appendix 
was solved with k c = 0. With this change, all the results 
obtained in appendix A of reference 8 are applicable here. 
The stability criterion for combined compression and shear 
is therefore the same in form as that for shear alone, given 
by equation (A-19) of reference 8, which is 

2a0^4y 3 — ^ (cosh 2 a cos 2j8— cos 4y)— a 2 ) 

— (/3M-0 2 )*— (4ry* — 0 l +a J ) sinh 2a sin 20 

+e[a(4rf+# > +0 > ) cosh 2 a sin 20 

+0(4y* — a* — 0 s ) sinh 2 a cos 20— 4a0y sin 4y]=0 (B7 

The relation between k, and a, 0, and y is also the same in 
form as that in equation (A-23) of reference 8, namely, 


7. . 87(0*+ 0 1 ) 

, xb 
®X 


(B8) 


In the present report, however, a and 0 have the following 
values: 


where 


V 0+ V w+c) '-re(T)‘ 

^=V-»+V (7!+<!)! -n(v)‘l 

, , 1 / x&Y x 4 , 
c- ^ + 4 (.X/ 8 kc 


(B9) 
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As in reference 8, the restraint coefficient e is defined herein 
by the relation 

Sb 

e= V "" 

where S is the ratio of a sinusoidally applied moment to the 
resulting sinusoidally distributed rotation of the restraining 
element measured in radians. 

Evaluation of k, corresponding to a selected value of k e . — 
The procedure for evaluating k„ after values of k e and e 
have been chosen, is as follows: A value of 6/X is selected; a 
series of values of y are assumed until one is found that, 
together with the corresponding values of a and j8 computed 
from equation (B9), satisfies equation (B7); k, is then com- 
puted from equation (B8), Another value of 6/X is selected; 
a new set of values of y, a, and 0 is found that satisfies the 
stability criterion; and a new value of k, is computed. The 
entire process is repeated until the minimum value of k, 
can be found from a plot of k, against b/\. When e is a 
function of b/\ e must be reevaluated each time a different - 
value of 6/X is selected. The minimum value of k, and the 
chosen value of k c , when inserted in equations (B2), give a 
critical combination of shear and direct stress. 

Evaluation of k, when k c has value corresponding to 
buckling as Euler strip. — One critical combination of shear 
and compression is simply &,=0 and k c equals the value 
corresponding to buckling as an Euler strip. The curves 
giving critical stress combinations, however, did not appear 
to be approaching this point as their construction pro- 
gressed. It was therefore necessary to determine whether 
values of k, other than zero are critical when the Euler com- 
pressive stress is reached. The determination of k, when 
k c reaches the value at which the plate buckles as an Euler 
strip requires special treatment, because k„ given by equa- 
tion (B8), becomes indeterminate when the wave length be- 
comes infinite as suggested by the energy solution. The 


result that X becomes infinite when k c takes its Euler value 
is readily checked from equation (Cll) for the special case 

of e=0; for this case p=g=r=^ and (£ c )<w=l. From 

equation (B8) it is clear that, if k, is to remain finite when 
the wave length approaches infinity, either y— >0, case (a), or 
a 2 — f j9 2 — >0, case (b). 

For case (a), when «=0, it follows from equations (B9) 
that, to small quantities of the second order. 


where 




(BIO) 


If the values of a and (3 from equations (BlO) are substi- 
tuted in equation (B7) and the resulting equation is ex- 
panded, with only the lowest powers of u and y retained, 


„ 4tt*U 2 

64— 6tt 2 


(Bll) 


Substitution of the values of a and jS from equations (BlO) 
and (Bll) in equation (B8) gives as the final result for «=0 

7. — ± -2t_ .. _ 

* ± V64 - 67T* 

Case (b) can be analyzed by a similar method, but the 
analysis is quite complicated because terras of third order 
must be retained. For e=0 and t— °o , case (a) and case 
(b) were found to lead to exactly the same result for k,. 
A value of k, other than zero when k c takes its Euler value 
may be found in the same manner for other values of edge 
restraint.. For any value of the restraint coefficient e, 



kc * 

(erV^+2, rVF e + WE> 

) sin vV^c+ 


2 ~~ 

\ COS 7T (COS v^k c — 1) 


1 


isin x-y 

r,+( 

3 **,. 1 

. 4 e 1 6 4 ) 

|COS rfVfcc+l 

(4+8t+2e ! + 

2^)(cosWF-1) 


(B12) 


where k c has the value corresponding to Euler strip buckling, 
at this value of e. The relationship between this Euler value 
of k c and the corresponding e is given by the equation 
(from reference 9) 

— _ 

tan^ - • ' ~ 

When k c reaches the Euler value, the critical shear-stress 
coefficient k, can therefore be either 0 or the value given by 
equation (B12). The conclusion that k, can also have any 
value between these limits is plausible on the basis of the 
following physical considerations. The shear stress does no 


work during buckling when the stress condition is such that 
the plate buckles with an infinite wave length. The effect 
of shear, furthermore, is to reduce the wave length to a value 
of the order of the width of the plate. The wave length at 
the time buckling occurs is infinite, however, when the plate 
is either in pure compression or at the value of k, satisfying 
equation (B12). This fact means that, for values of k, 
between 0 and that given by equation (B12), the shear 
stress is not great enough to force buckling in short waves 
and therefore does not assist in producing buckles. In this 
range of shear stress the compressive stress necessary to 
produce buckling is, consequently, the Euler stress. 



APPENDIX C 


SOLUTION BY ENERGY METHOD FOR EDGE RESTRAINT INDEPENDENT OF WAVE LENGTH 


The critical stress is determined on the basis of the prin- 
ciple that the elastic-strain energy stored in a structure 
during buckling is equal to the work done by the applied 
loads during buckling. If the structure under consideration 
is an infinitely long plate under combined shear and edge 
compression with edges elastically restrained against rota- 
tion, this equality may be written 


Tc+T.^+V, (Cl) 

In reference 1 an energy solution was given for the type 
of loading shown in figure 1 (b). The deflection function 
used in reference 1 is also suitable for application to the 
solution of the type of loading shown in figure 1 (a), which 
is the loading considered in the present paper. The values 
for T„ Vi, and Vi may accordingly be taken directly from 
reference 1, but T e must be recomputed to apply to the case 
of transverse compressive stress. 

The following substitutions are used to transform the 
energy expressions from the oblique coordinates of reference 
1, the left-hand terms of the equations, to the rectangular 
coordinates (fig. 6) used in the present paper, the right-hand 
terms of the equations: 

V 

bTb\ (02) 

bi cos <t>—b ) 

For brevity the following notation is also adopted 

tan 4>=d 


1 

/ 

_J_J Ul 

V / 

y/ ^ 

b / 

? / L_ 

■ 7 *=^ 

i 

i 

V A 

t — r^~i — r — 

“f — 1 — t — I — 1 — t 


Floral 6.— Coordinate system used and wave form assumed for energy solution In appendix C ; 
Inclined lines Indicate nodal lines of buckles. 


By use of equations (C2), the expressions from reference 1 
that are used in the present paper may be rewritten as 
follows : 


T,=w<? 


N xl/ x 2 bd j. 
2 \ Jl 


(C3) 


[0\l + Wi+2(1 + 3^)/s+Q) 2 ft] (04) 


F,=«b* 


TT*D\t 

~w~ 


(05) 


where 

fl= (S + §~ J 5 ) '+(*-?) e4 i 

/s= (^~^) e2+ G~^) e+ i 

and £ is the restraint coefficient defined in appendix B. 

The work done by the compressive force per half wave 
length may be written 

T --k N -fSXW ixiy (C6) 

2 2 


As in reference 1, the assumed deflection function is 

W=W 0 Jjjr (Js“£)+( 1 +§) cos x] cos X ^ ( C7 > 

When the expression for w from equation (07) is substi- 
tuted in equation (C6) and the indicated operations are per- 
formed, 

Tc=Wo*^(jJ*+Vi) ( 08 ) 

Values from equations (C3), (C4), (C5), and (C8) are now 
substituted in the buckling equation (Cl). The use of the 
equations 

]y _j t 

iy xv — ic, 


N u =k, 


b 2 

t^D 

V 


eliminates N xv and N„. The resulting equation gives the 
critical combination of stresses and may be written as 

*‘ = 2^[JS ( 1 + Wi+F/*+ 2 O + 30 2 )/* 

+Pp-^(p/ a +^)] (09) 

This equation shows that, for a selected value of k C) the criti- 
cal shear stress depends upon the wave length and the angle 
of the buckle. Since a structure buckles at the lowest stress 
at which instability can occur, k, is minimized with respect 
to wave length and angle of buckle. The minimum value of 
k, with respect to value of wave length is determined from 
the condition 

m-° <° 10 > 
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which gives (when e does not depend on wave length) 



d+e 2 )^ 


V' 


7T 


(Cll) 


Substitution of this value of wave length in equation (C9) 
gives 


k,=1 w i* [@+/3)-^]i i/2 +4 [2(i +3^- Wl ] 

(C12) 


The minimum value of k, with respect to angle of buckle is 
found from the condition 


be u 


which gives 



/i0rh/ s ) 

\—kJ% 

1/2 

+/* 



4 1 

-3 

'A 

W +3/,-hp 


(C13) 


(C14) 


If this value of 6 is substituted in equation (Cl2), the final 
result is the following interaction formula, in which k, is 
given in terms of k e and the edge restraint e: 


*.*■=401+301 (Ct-kc) + (4 C 2 -k c ) V2(2Ci - C 2 k c ) (Cl5) 


where 


and 
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TABLE I.— VALUES OF ft, AND R c WITH CORRESPONDING 
COMPUTED VALUES OF ft, AND R, 


k. 

k. 


R, 

Exact 

Energy 

Rt 

>Vr 

Exact 

Jik' 

i^O 

1.00 

r o 
\ ase 

0 

8.16 

} 1.00 

f 0 

X .'335 

.99 

3.09 

3.39 

.99 

.578 

.95 

3.36 

8.67 

.95 

.629 

.90 

8.67 

3.88 

.90 

.609 

.80 

a 87 

4.20 

.80 

.726 

.60 

4.35 

4 66 

.60 

.814 

.40 

4.70 

6.04 

.40 

.881 

.30 

4.88 

5.20 

.30 

.913 

.20 

6.03 

5.36 

.20 

.942 

.10 

a 19 

5.51 

.10 

.972 

0 

a 34 

aoe 

0 

-1.000 

-1.00 

6.56 

6.89 

-1. 00 

1. 229 

-3.00 

7.63 

7.91 

-200 

1.428 ’ 

-5.00 

10.22 

10.47 

-400 

1.914 

«-2 

L 67 

X 2.77 

0 

2.06 

} 1.00 

/ 0 
X .457 

1.33 

4.07 

4.82 

.80 

.872 

.67 

6.20 

5.44 

.40 

.868 

0 

6.06 

6.30 

0 

1.000 

-8.34 

9.28 

9.60 

-2.00 

1.627 

-8.34 

12.95 

ia 16 

-400 

2137 

«»6 

2.30 

/ 0 
l 2.76 

0 

ao3 

} 1.00 

f 0 

X .413 

L84 

4.27 

4.43 

.80 

.639 

.92 

6.62 

a8i 

.40 

.841 

0 

6.68 

ass 

0 

LOGO 

-4.59 

10.67 

10.88 

-200 

1.597 

-11. 49 

ia28 

iaR2 

-400 

2237 

€-10 

2.86 

1 0 
X 178 

0 

8.18 

} 1.00 

f 0 

X .382 

2.28 

4 49 

460 

.80 

.618 

1. 14 

6.04 

6.19 

.40 

.831 

0 

7.27 

7.43 

0 

1.000 

-5.71 

11.87 

12.13 

-2.00 

1.C33 

-14.27 

17.29 

17.60 

-5.00 

2 378 

(■ CO 

4.00 

/ 0 
1 3.27 

0 

3.81 

} 1.00 

f 0 

X .364 

aso 

4.27 

4 34 

.96 

.476 

aeo 

4.73 

476 

.90 

.627 

a20 

6.42 

445 

.80 

.004 

2.40 

a52 

457 

.60 

.726 

L 60 

7.42 

7.52 

.40 

.826 

.80 

8.24 

437 

.20 

.917 

0 

a98 

9.15 

0 

1. 000 

-4.00 

12.21 

12.62 

-L00 

1.300 

-a oo 

14 95 

14 39 

-2 00 

1.665 

-20.00 

21.90 

22,67 

-400 

2439 




